LOW MACH NUMBER FLOWS, AND COMBUSTION 



THOMAS ALAZARD 

Abstract. We prove uniform existence results for the full Navier-Stokes 
equations for time intervals which are independent of the Mach num- 
ber, the Reynolds number and the Peclet number. We consider general 
equations of state and we give an application for the low Mach number 
limit combustion problem introduced by Majda in |18| . 



1. Introduction 

For a fluid with density velocity u , pressure P, temperature T, internal 
energy e, Lame coefficients C, r] and coefficient of thermal conductivity /c, 
the full Navier-Stokes equations, written in a non-dimensional way, are 

' dtp + div(pu) = 0, 

VP 

(1.1) I dt{pv) + di\{pv ®v) + — = ^i{2di\{C,Dv) +V{r]d\Yv)), 
^ dt{pe) + div(pwe) + P diw = Kdiv(A:VT) + Q, 

where e G (0, 1], (/i, k) € [0, 1]^ and Q is a given source term (see PITIITHIITS] ). 
In order to be closed, the system is supplemented with two equations of state, 
so that p,P,e,T are completely determined by only two of these variables. 
Also, it is assumed that (, rj and k are smooth functions of the temperature. 

This paper is devoted to the asymptotic limit where the Mach number e 
tends to 0. We are interested in proving results independent of the Reynolds 
number 1/p and the Peclet number 1/k. Our main result asserts that the 
classical solutions of ()1.1() exist and are uniformly bounded on a time interval 
independent of e, ^ and k. 

This is a continuation of our previous work [1^ where the study was re- 
stricted to perfect gases and small source terms Q of size 0(e). We refer 
to the introduction of ^ for references and a short historical survey of the 
background of these problems (see also the survey papers of Danchin j^l, 
Desjardins and Lin Gallagher [13], Schochet [M] and Villani ^). 

The case of perfect gases is interesting in its own: first, perfect gases 
are widely studied in the physical literature; and second, it contains the 
important analysis of the singular terms. Yet, modeling real gases requires 
general equations of state (see |3 Ej)- Moreover, we shall see that it is 
interesting to consider large source terms Q for it allows us to answer a 
question addressed by Majda in ^H] concerning the combustion equations. 
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1.1. The equations. To be more precise, we begin by rewriting the equa- 
tions under the form L(n, dt, dx)u + s~^S{u, dx)u = 0, which is the classical 
framework of a singular limit problem. 

Before we proceed, three observations are in order. Firstly, for the low 
Mach number limit problem, the point is not so much to use the conservative 
form of the equations, but instead to balance the acoustics components. This 
is one reason it is interesting to work with the unknowns P,v, T (see jl8j). 
Secondly, the general case must allow for large density and temperature 
variations as well as very large acceleration of order of the inverse of the 
Mach number (see Section 5 in Since dtv is of order of e~^VP, this 

suggests that we seek P under the form P = Cte + 0(e). As in [20j, since 
P and T are positive functions, it is pleasant to set 

(1.2) P = Pe^P, T = ^e^ 

where P_ and T are given positive constants, say the reference states at 
spatial infinity. Finally, the details of the following computations are given 
in the Appendix. 

From now on, the unknown is (p, v, 0) with values in M X M'^ X M. We are 
interested in the general case where p and 9 are uniformly bounded in e (so 
that VT = 0(1) and dtv = 0{e-^)). 

By assuming that p and e are given smooth functions of (P, T) , it is found 
that, for smooth solutions of ()1.1|) . {P,v,T) satisfies a system of the form: 

( a{dtP + V ■ VP) +dwv = k(3 dw{kVT) + pQ, 



(1.3) 



VP 

p{dtv + V ■ Vv) + = fi{2dw{CDv) + V(r/divi;)), 
t j{dtT + V ■ VT) + div f = k5 div(A;Vr) + 5Q, 



where the coefficients a, /?, 7 and 5 are smooth functions of {P,T). Then, 
by writing dt^xP = ^Pdt,xP, dt^xT = Tdt,xO and redefining the functions 
C, and T], it is found that {p,v,6) satisfies a system of the form: 



(1.4) 



9i{^){dtP + v-Vp) + ^divv= '^xi{(t>) div{k{9)Ve) + ^Xi{(t>)Q, 
92(.4>){dtv + V ■ Vv) + ^Vp = pB2{4>, dx)v, 



^ gsmdtO + vve) + dwv = KX3{<P) dw{kie)V9) + x3i<l>)Q, 

where := {9,ep) and 52(0,9.) = X2{<P) dw{Ci9)D-) + X2{(l>)V{rj{9) div ■). 

We are now in position to explain the main differences between ideal gases 
and general gases. Firstly, we note that the source term Q introduces an 
arbitrary unsigned large term of order of 1/e in the equations. Secondly, to 
emphasize the role of the thermodynamics, we suppose now that Q = and 
we mention that, for perfect gases, the coefficient Xi{4') is a function of ep 
alone (see Prop osit ion I A . 8|) . Hence, for perfect gases, the limit constraint is 
linear in the sense that it reads divide = with Ve = v — KXi{0)k{9)V9. By 
contrast, for general equations of state, the limit constraint is nonlinear. 
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1.2. Assumptions. To avoid confusion, we denote by (■!?, p) E M'^ the place 
holder of the unknown (9,ep). Hereafter, it is assumed that: 

(HI) The functions ij and k are C°° functions of G M, satisfying A; > 0, 

C > and r/ + 2C > 0. 
(H2) The functions gi and Xi = 1)2,3) are positive functions of 
p) gM?. Moreover, 

Xi < X3, 

and there exist two functions F and G such that {'d, p) i-^ {F{-&, p), p) 
and ("i?, p) ^ {'&,G{'&, p)) are C°° diffeomorphisms from M? onto M^, 
F(0,0) = G(0,0) = and 

OF dF dC dG 

Remark 1.1. Assumption (lEl) is used to prove various energy estimates. 
The main hypothesis is the inequality xi < Xs- Appendix IX} it is proved 
that the inequality xi < Xs holds whenever the density p and the energy e 
are C°° functions of (P, T) G (0,+oo)^, such that p > and 

^\ P^^T^= ^— ^ n dp ^ de dp 

^ ' dP dT ^ dp' dP ' ar ' drop dP or' 

1.3. Main result. We are interested in the case without smallness assump- 
tion: namely, we consider general initial data, general equations of state and 
large source terms Q. To get around the above mentioned nonlinear features 
of the penalization operator, we establish a few new qualitative properties. 
These properties are enclosed in various uniform stability results, which as- 
sert that the classical solutions of (|1.4j) exist and they are uniformly bounded 
for a time independent of e, p and k. We concentrate below on the whole 
space problem or the periodic case and we work in the Sobolev spaces H'^ 
equipped with the norms 11^11^:^^ := ||(/ — A)°"/^n||^2- 

The following result is the core of all our other uniform stability results. 
On the technical side, it contains the idea that one can prove uniform esti- 
mates without uniform control of the norm of the velocity v. 

Theorem 1.2. Let d = I or d ^ 3 and N 9 s > 1 -|- d/2. For all source 
term Q = Q{t,x) £ C^{R x M'^) and all Mq > 0, there exist T > and 
M > such that, for all {£,p, n) E (0, 1] x [0, 1] x [0, 1] and all initial data 
{po,vo,eo) e H'+\R'^) satisfying 

(1.6) IKVpo, Vfo)||j^.-i + 11(^0, epo, £1^0)11^^+1 ^ Mq, 

the Cauchy problem for ()1.4() has a unique classical solution (p, v, 9) in 
G^{[Q,T];H'+^{W^)) such that 

(1.7) sup \\{Vp{t),Vv{t))\\H.-, + \\{9{t),ep{t),ev{t))\\Hs ^ M. 
A refined statement is proved in Section |21 

A notable corollary of Theorem ll.2l is Theorem 14. 11 which is the requested 
result for application to the low Mach number limit. Detailed discussions of 
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the periodic case and the combustion equations are included in Sections El 
and El The assumption d 7^ 2 is explained in Remark 12.61 

2. Preliminaries 

In order not to interrupt the proofs later on, we collect here some esti- 
mates. The main result of this section is ProDosition l2.41 which complements 
the Friedrichs-type estimate 

(2.1) llViill^s ^ lldivfllji^s + ||curlt;||j:^s , 

which is immediate using Fourier transform. We prove a variant where div v 
is replaced by div{pv) where p is a positive weight. 

Notation. The symbol < stands for ^ up to a positive, multiplicative 
constant, which depends only on parameters that are considered fixed. 

2.1. Nonlinear estimates. Throughout the paper, we will make intensive 
and often implicit uses of the following estimates. 

For all fj ^ 0, there exists K such that, for all u,v £ D H'^{R''-), 

(2.2) ll'ut^llH'' ^ ^ II ^IIl^ II^II^ct ~t~ K \\u\\jja ll^ll/^oo • 

For all s > d/2, fii ^ 0, (T2 ^ such that ai + a2 ^ 2s, there exists a 
constant K such that, for all u G F'*"'^i(M'^) and v G H^~''^{R'^), 

(2.3) 1 1 1 1 ^"'''1 ^ ^ W'^W JJ^~^2 ' 

For all s > d/2 and for all C°° function F vanishing at the origin, there 
exists a smooth function Cp such that, for all u G H''{W^), 

(2.4) \\F{u)\\hs^Cf{\\u\\l^)\\u\\hs. 

2.2. Estimates in W^. Consider the Fourier multiplier VA~^ with symbol 
— i^/lCl^- This operator is, at least formally, a right inverse for the diver- 
gence operator. The only think we will use below is that VA~^n is well 
defined whenever u = U1U2 with ui,U2 G L°° n H'^{W^) for some o" ^ 0. 

Proposition 2.1. Given d ^ 3 and cr G M, the Fourier multiplier VA^-*^ is 
well defined on (W^) D H'^ (W^) with values in i^'^+^(M'^). Moreover, there 
exists a constant K such that, for all u G L^{W^) n H^'i/W^), 

(2.5) ||VA~^ti||^,+i ^ K + K WuWjja . 

Proof. Set (0 := (1 + 1^1^)^^^- It suffices to check that the L^-norm of 
{{iY^^I is estimated by the right-hand side of ^T^. To do that 

we write 

/ ^7^|S(0l''ie<l|n||i., / AiL^|2(^)|2d^<||^||2^^, 

where we used 1/ G L^{{\i\ ^ 1}) for all d ^ 3. □ 

The next proposition is well known. Its corollary is a special case of a 
general estimate established in j^. 
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Proposition 2.2. Given d ^ 3 and s > d/2, there exists a constant K such 
that, for all u G H'{W^), 

(2.6) \\u\\j^oo ^ K\\Vu\\^s-i ■ 

Proof. Since if (M"') ^ L'^iW^), it suffices to prove the result for u in the 
Schwartz class 5(M'^). Now, starting from the Fourier inversion theorem, 
the Cauchy-Schwarz inequality yields the desired estimate: 

11-11.^ ^ (/ ^^^^(iru)''' (/ (0^^-^^ dij'' < ||Vn||,._. . 

Corollary 2.3. Given d ^ 3 and N 9 s > d/2, there exists a constant K 
such that, for all ui,U2 € H^{R.'^), 

(2.7) II ^i^llVnill 

Proof. One has to estimate the L^-norm of d^{uiU2), where a G N'^ satisfies 
|a| ^ s. Rewrite this term as uid^U2 + [d°,ui]u2. Since the commutator 
is a sum of terms of the form dxUid2u2 with (3 > 0, the product rule H2.3() 
implies that 

(2.8) ll[5:,ni]n2||i2<||Vni|| 

Hs-i ll'"2||/fs . 

Moving to the estimate of the first term, we write 

\\uid^U2\\L2 ^ IKlllioo \\u2\\h= < ||Vui||^s-i ||^i2|k=. □ 

2.3. A Friedrichs' Lemma. With these preliminaries established, we are 
prepared to prove the following: 

Proposition 2.4. Let d ^ 3 and N 9 s > d/2. There exists a function C 
such that, for all cp G ii^+i(R'^) and all vector field v G H''+^{R'^), 

(2.9) l|V?;||^. ^ C ||div(e'^w)||^. + C ||curl7;||^. , 
where C := (1 + ||(/?||^,+i)C(||v9||j:^, , || V(/?||^oo) • 

Proof. For this proof, we use the notation 

R = ||div(e'^t;)||^. + ||curlt;||j:^. , 

and we denote by various constants depending only on ||</7||j:^s + ||V(^||^oo. 

All the computations given below are meaningful since it is sufficient to 
prove ()2.9|) for C°° functions with compact supports. We begin by setting 

V = V + \7A'^{\7(p ■ v). 

The reason to introduce v is that 

e''' diviJ = div(e'''t'), curl-u = curlu . 

Hence, by using 1)2. we have 

(2.10) llVi^ll^. ^ ||e~'^div(e^w)||^, + ||curl7;||^. ^ C^R. 



6 



T. ALAZARD 



The proof of (|2.9j) thus reduces to estimating vi := v — v, which satisfies 

div(e'^vi) = —e^Vif ■ v, curlui = 0. 

Again, to estimate vi we introduce vi := vi + VA~^(Vip ■ vi), which solves 

divui = —Vif ■ V, curltJi = 0. 

The estimate ()2.ip imphes that HViJillji^s ^ ||V(/7 • By using H2.1U() 

and the product rule H2.7|) . applied with ui =v and U2 = V(/3, we find that 

(2.11) llViJill^, ^ \\V(p-v\\jjs < \\(p\\H^+i ||ViJ||j^,_i ^ \\ip\\j^s+i C^R. 
Hence, it remains only to estimate V2 = vi — vi, which satisfies 

div(e'^W2) = —e'^Vip • vi and curlu2 = 0. 
To estimate V2 the key point is the estimate 

(2.12) \\V^ ■ viW^d* i^C^R, with d* = 2d/{d + 2). 

Let us assume ()2.12() for a moment and continue the proof. 

The constraint curlt>2 = implies that V2 = V^*, for some ^ satisfying 

div(e'^V^) = -e'^Vip-vi. 

This allows us to estimate V^' by a duality argument. We denote by ( • , • ) 
the scalar product in and write 

Denote by d the conjugate exponent of d*, d = d* /{d* — 1) = 2d/{d — 2). 
The Holder's inequality yields 

The first factor is estimated by means of the claim (|2.12|) . In light of the 
Sobolev's inequality ||^||^d ^ I|V^IIl2) obtain 

By using the elementary estimate ||V^||^2 ^ ||e~'''|lLoo (e'^V'I', V^f), we get 

(2.13) \\v2\\l2 = \\V^\l2 ^ C^R. 
The end of the proof is straightforward. We write 

= e"^ div(e^V^') - Vy? • V^' = -Vip -vi-VLp- V^, 

to obtain, for all o" G [0, s — 1], 

l|V^'||^^.+i < ||V*||^2 + ||A^'||^. < llVvp-^Jill^. + (l + ||v;||^.)l|V^'||^^. . 

To estimate the first term on the right-hand side, we verify that the analysis 
establishing (|2.7|) also yields 

||V(/?-^i||^.-i < W^pWhs ||Vui||^s-i ^ C^R, 

hence, by induction on a, 

l|V*||^. ^C^R + C^\\V^\\l2. 
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Exactly as above, one has 

llV^'lljy.+l < ||V^'||i2 + IIA^-II^, < ||V(/^-?^l||^. + (l + ||(/j||j:^.+ l)||V^'||^. 

As a consequence, we end up with 

||V^'||j:^. + l ^ M\hs+1 (C<^-R + C<p||V^||^2). 

Therefore, the estimate ()2.13() imphes that 

||'y2||_H-s+i = llV^fllji^s+i ^ \\ip\\jjs+i C^R. 



By combining this estimate with (|2.11|) . we find that 

From the definition of vi and (|2.1()|) . we obtain the desired bound (|2.9j) . 
We now have to estabhsh the claim (|2.12|) . 

With d = 2d/{d — 2) as above, the Sobolev's inequality and (|2.1Up imply 
that 

ll^^llid < l|V^J|li2 ^ C^R. 
On the other hand, the Holder's inequality yields 

^ l|V(/?IL2 ll^ll^d , with 5 = = 
By interpolating this estimate with ||V(^ • v\\^-^ < ||V99||j;^oo we obtain 

Because curl?;i = 0, one can write vi = V^'i for some function ^'i sat- 
isfying A\I'i = — V(/7 • V. Hence, the Calderon-Zygmund inequality and the 
previous estimate imply that 

||V7;i||^. = p^^iW^, < WA^iWl^ ^ C^R- 
Therefore, the Sobolev's inequality yields 

Sd d 

HW.n^C^R, Withi?=^ = ^, 
hence, exactly as above, the Holder's inequality gives 

2D 2d 



(2.14) yp£[d,d], \\e^Vip-vi\\^,^C^R, with d : _ ^ 

The key estimate (|2.12|) is now a consequence of the previous one. Indeed, 
the estimate (|2.14p applies with p = d* = 2d/{d + 2) since 

, 2d 2d 2d - 

yd ^3, d= — ^ ^ = d. 

- 3d-A d + 2 d-2 

This completes the proof of 1)2. 9() . □ 
For later references, we will need the following version of (|2.9|) . 
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Corollary 2.5. Let d = 1 or d ^ 3 and'N B s > d/2. There exists a function 
C such that, for all G and all vector field v G H''+'^{W^), 

(2.15) ||Vt;||^. ^ CiMHs+i){\\dWv\\Hs + ||curl(e'^t;)||^.). 

Proof. The case d = 1 is obvious. If d ^ 3, Proposition 12.41 (applied with 
{(p,v) replaced with {—ip,e'^v)) yields 

l|V(e'^t;)||^. ^ Cdlv^ll^.+OdldivHI^. + ||curl(e'^t;)||^.). 

Hence, to prove (|2.15jl we need only prove that 

(2.16) ||V7;||^. ^ C(||<^||^.+i) ||V(e^7;)||^. . 

To do that we write diV = e~^di{e^v) — {e~^diLp){e^v). The usual product 
rule (|2.3|) implies that the norm of the first term is estimated by the right- 
hand side of 1)2. 16(1 . Moving to the second term, we use the product rule 
(EH) to obtain \\{e~'Pdi^){e'fv)\\jjs < (1 + M\hs) \\dM\H^ ||V(e^i')||j;^._i . 
This proves the desired bound 1)2. 16() . □ 

Remark 2.6. The fact that Theorem 11.21 precludes the case d = 2 is a 
consequence of the fact that we do not know if (|2.15|) holds for d = 2. 

3. Uniform stability 

In this section, we prove Theorem 11.21 We follow closely the approach 
given in we recall the scheme of the analysis and indicate the points at 
which the argument must be adapted. 

Hereafter, we use the notations 

a := (e, fi,K) £ A := (0, 1] x [0, 1] x [0, 1], u := ^^uTk, 
||u||^<T+i := ||ti||j|^<T + a ||n||j|^CT+i (a ^ 0, a £ M). 

Step 1: a refined statement. We first give our main result a refined form 
where the solutions satisfy the same estimates as the initial data do. Also, 
to prove estimates independent of fi and k, an important point is to seek 
the solutions in spaces which take into account an extra damping effect for 
the penalized terms. 

Definition 3.1. Let T > 0, a = {e, fi, k) G [0, 1]^ and set v = yj ii + k. The 

space X^{T) consists of these {p,v,9) G C°([0, T]; i7^(M'^)) such that 

u{p,v,d) G C\[0,T];H'+\R'^)), {fiv,Ke) G L\0,T; H'+\R^)). 
The space A'^(T) is given the norm 

II(p>''^>^)IIa'^(T) ■= ll(VP'Vw)||^op(j^.-i) + \\{e,ep,ev)\\^^^^s+i^ 

+ VJ^ \\^'^\\lI{hi+') + ^ \\^^\\Lum+') 

+ VJ^^\Np\\lI{h^) + ^ l|divw||^2 (j:^.) , 
with (X) denoting the norm in LP(0,T; X). 
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The hybrid norm was already used by Danchin in 

For the study of nonhnear problems, it is important to relax the assump- 
tion that Q £C§°. 

Definition 3.2. The space -F" consists of these function Q such that, for 
aUn3m^s, d^Q £ C^{R; H'+^-^'^'{R'^)), where stands for C° n L°°. 

Given a normed space X, we set B{X; M) = {x £ X : ||x|| ^ M}. 

Theorem 3.3. Assume that d = I or d ^ 3 and let 'N 3 s > 1 + d/2. 
Given Mq > and Q £ , there exist T > and M > such that, for all 
a = {e,iJ,,K) G A and all initial data {po,Vo,d(}) G H'^~^^{M.'^) satisfying 



the Cauchy problem for H1.4() has a unique solution {p,v,0) G B{X^(T); M). 

This theorem implies Theorem 11.21 

Remark 3.4. A close inspection of the proof indicates that Theorem 13.31 
remains valid with p.lf) replaced by 

\\{Po,vo,6o)\\x^f{o) ■= ll(Vpo,Vt;o)||/^s-i + ||(6'o, epo, et^o)||/^j+i ^ Mq. 

Step 2: local well posedness. We explain here how to reduce matters 
to proving uniform bounds. To do so, our first task is to establish the local 
well posedness of the Cauchy problem for fixed a = {e, fi, n) G A. 

Lemma 3.5. Let d ^ 1, s > 1 + d/2 and a £ A. For all initial data Uq = 
{Po,vo,9q) G i7^(M'^), there exists a positive time T such that the Cauchy 
problem for (|1.4|) has a unique solution U = {p,v,6) G C^{[0,T]; H'^) such 
that U{0) = Uq. Moreover, the interval [0,T*), withT* < +oo, is a maximal 
interval of H"" existence if and only lim supj^^-* 11^^(^)11 \Yi,oc md\ — -|-CX). 

Lemma 13.51 is a special case of Proposition 14.51 established below. 
As in ^ 1201 , on account of the previous local existence result for fixed 
a £ A, Theorem II. 21 is a consequence of the following uniform estimates: 

Proposition 3.6. Lei d = 1 or d ^ 3, N 9 s > 1 + d/2 and Mq > 0. 

There exist a constant Cq and a non-negative function C(-) such that, for 
allT G (0,1] and alia G A, if {p,v,e) G C°°{[0,T];H'^{R'^)) is a solution 
of (|1.4|) with initial data satisfying (|3.H) . then the norm 0,a{T) := {{UW^s/j'^^ 



Notation 3.7. From now on, we consider an integer s > 1 + d/2, a fixed 
time < T ^ 1, a fixed triple of parameters a = (e, fi, k) G ^, a bound Mq, 
a fixed smooth solution U = {p,v,9) G C~([0, T]; of Jlil) with 

initial data satisfying (|3.1j) and we set 



(3.1) 



\\{Vpq,Vvq)\\hs-i + \\{9Q,epQ,evQ)\\Hs+ 




naiT) ^ CQexp{iVT + e)C{naiT))). 



n := \\u 



XS{T) ■ 
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With these notations, Prop osition ll-? . 61 can be formulated concisely as follows: 
if d 7^ 2, there exist constants Co depending only on Mq and C depending 
only on Vl such that 

Hereafter, we use the notations (j) := {9,£p) and v := 

Notation 3.8. For later application to the nonlinear case when Q = F{Y) 
for some unknown function Y, we also give precise estimates in terms of 
norms of Q. For our purposes, the requested norm is the following: 

(3.3) S:= V \\{I-{eufAr'''^'{e{dt + vV)y"'Q 



L°°(0,T;_f/J+i-'") 



L°°(0,T;//J+l-'") 



Remark 3.9. To use nonlinear estimates, it is easier to work in Banach 
algebras. If d ^ 3, Proposition 12.21 shows that we can supplement the 
estimates with L°° estimates for the velocity: it suffices to prove ()3.2() with 
C{Qa{T)) replaced by C(S7+(r)) where 0+(r) := ^^(r) + ||w||^cx,((o,t)xIR'*)- 
Similarly, if d ^ 3, all the estimates involving the source term Q remain 
valid with T, replaced by 

||(/-M'A)-™/2(e5trQ 

Step 3: An energy estimate for linearized equations. A key step in 
the analysis is to estimate the solution {p,v,6) of linearized equations. As 
will be apparent in a moment, a notable fact is that we can see unsigned large 
terms e~^f^{t, x) in the equations for p and v as source terms provided that: 
1) they do not convey fast oscillations in time: dtf^ = 0{1); 2) it does not 
implies a loss of derivatives. To be more precise: in the nonlinear estimates, 
we will see the term £~^Xi{4')Q ^ source term. Similarly, we can see terms 
of the form £-^F{£p,e, ^J1^\J9) as source terms. As a result, it is sufficient 
to consider the following linearized system: 
(3.4) 

Ik ~ 
gi{(t)){dtp + V ■ Vp) + -diviJ- -div(A:i((/>)V6') = Fi, 

g2{4'){dtv + vVv) + ^Vp - fiB2{(l), d^)v = F2, 
^ gm{dtO + v-Ve) + G{ct>, Vct))-v + divv- Kx^icj^) div(A:(0)V^) = F3, 

where the unknown {p,v,9) is a smooth function of {t,x) G [0, T] x M^. 
The following result establishes estimates on 

\\iP,V,9)\\a,T ■■= ll(P,^')llL-(iTiJ + \\G\\l^(HI) 

(3-5) + V^\\V9\\l2^^hI) + VJ^ miLliH^^) 

+ ^/JT^ II Vp||i2 (^2) + ^ lldivl^ll ^2^(^2) , 

in terms of the norm ||(p, i;, 0)||a^o ||(P;''^)(0)lli^i +11^(0)11^1 °^ data. 
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Theorem 3.10. Let d ^ 1 and assume that G, ki and are C°° functions 
such that, for all p) e M^ < ki{i}, p) < xsi^, p)H^)- Set 

Ro ■■= UiO)\\Hs-^ ^ R-= sup \\{(k,dt<f) + v-V(t),V<f),iyV^(k,Vv){t)\\ 

iG[0,T] 

There exist constants Cq depending only on Rq and C depending only on R 
such that, 

\\{j>,vMa,T^Coe^''\\iPo.vM\\afi + C || (Fi , F2 ) || +11^^311^1 dt. 

Jo 

In we established the previous theorem with Rq and R replaced by 
Ro = Um\L^, R'= sup \\{^,dtcl),v,V^,uV^^,Vv){t)\\^^. 

t€[0,T] 

To prove the above variant, we need only check two facts. Firstly, in the 
proof of Theorem 4.3 in the terms dt(p and v always come together within 
terms involving the convective derivative + v ■ Vcf). 

Secondly, we have to verify that the norms of the coefficients ((7, ((/>),...) 
are estimated by constants of the form Cqc"^^ . In ^ we used the estimate 

sup ||F(<A(t))||ioo ^ ||F(</.(0))||^oo+r sup \\dtFm))\\^^^C'^ + TC\ 

iG[0,T] tG[0,T] 

for some constants depending only on R'q and R! . Here, based on an usual 
estimate for hyperbolic equations, we can prove a similar bound: 

Lemma 3.11. Let F G C°°(M2) he such that F(0) = 0. There exist con- 
stants Cq depending only on Rq and C depending only on R such that, for 
allt£ [0,r], ^ Coe^^. 

Proof. Since s — I > d/2, the Moser's estimates (|2.'2() and H'2.4() imply that 
there exists a function C depending only on the function F such that 

\\{dt+vv)FmHs-i 

^ (1 + ||F'(</.) - F'(0)||^,_0 11(9^ + V ■ V)m\Hs-^ , 

and \\Fim)\\H^-i^C{\\m\\H^-i)- 

Hence, the desired estimate follows from the following estimate: there 
exists a constant V depending only on || Vt>||^«)j|^s-i such that 

sup ^ \\F{m)\\Hs-i+TV sup||(9t + t;-V)F(0(t))||^._i. 

te[o,T] te[o,T] 

To prove this result we set S := d'^F{(f)) where a G N*^ is such that \a\ ^ s—1. 
Then u solves 

dtu + vVu = f := {{dt + V ■ V)F{<P)) + [v, 8^] ■ VF{<P). 
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Since s — 1 > d/2, the product rule ^l.'A^ implies that 
\\[v,d^]-VFmL^ < E \\d^vd2VF{<p)h2 

/3+7=a, /3>0 

^ X] ll^f^'llH^-i-d/si-i) ll^2V^lliy»-i-(i-Ti+i) , 

/3+7=a, /3>0 

hence, ||/||^. < 11(9* + t; • V)F(</.)||^._, + ||V7;||^.-i ||i^(<A)llH-i- 
We next use an integration by parts argument yielding 

j^\\n\\l2 ^{l + \\divv\\^^)\\u\\l, + \\f\\l,. 

The Gronwall's Lemma concludes the proof. □ 

Step 4: High frequency estimates. On the technical side, the estimate 
of the derivatives is divided into four steps. Most of the work concerns the 
separation of estimates into high and low frequency components, where the 
division occurs at frequencies of order of the inverse of ez^ where = ^ + k. 
We begin by estimating the high frequency component 



n 



HF 



:= \\{I-Je,)U\\^^^T), 



where {Jh\h €z [0, 1]} is a Friedrichs mollifiers: Jh = j{hDx) is the Fourier 
multiplier with symbol j{h^) where j is a function of ^ G M"^, satisfying 

O^J^l, Xe) = 1 for ICK 1, Xe) = for 1^1 ^ 2, = j {-^ . 

Proposition 3.12. Let d ^ 1. There exist constants Cq depending only on 
Mq and C depending only on 17, such that 



(3.6) ^ Coe^^ + VtC ||Q||^^(^.+i 



Proof. Introduce P := {I-Jeu)^' and U := {Pp,Pv,Pe). Then, U satisfies 
System (|^ with 

h{(t>):=xim{e), G(0,V0):=53( 
and F = {Fi,F2, F^f := /hf + /g + fx, where 



\, fx-=\ ' 



p{xmQ) ) \ 

and /hf is given by 

/i,HF= [gi{(t>),P]{dt + v-V)p + gi{c^)[v,P] -Vp -^[Bi{ct>,dx),P]e, 
/2,HF= [g2{(p),P]{dt + v-V)v + g2{(l))[v,P]-Vv -fi[B2{(l),dx),P]v, 
/3,HF = [93{(P),P] {dt + v- V)0 + g3{(l)){v; P}-Ve- A^[i?3(0, dx),P]e, 

where Bi{<j), d^) = xM) div(A:(6l)V-) {i = 1, 3), [A, B] = AB - BA and 
[v; P} ■V9:=v VPe + (Pv) • V9 - P{v • V^). 
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Estimate for /hf- We use the following analogue of Lemma 5.3 in there 
exists a constant K = K{d, s) such that 

II [/, P]u\\^, ^ ei^K WVfhoo \\u\\hs + ei^K ||V/||^. \\u\\l^ , 
II [/' PHhi ^ II^/IIl- MIh^ + vK ||V/||^. ||n||^^ . 

The fact that the right-hand side only involves V/ follows from the most 
simple of all the sharp commutator estimates established in jl7j : for all 
s > 1 + d/2 and all Fourier multiplier A{Dx) S OpSfg, there exists a 
constant K such that, for all / e H^iW^) and all u G i7*(M"'), 
(3.7) \\[f,A{D,)]u\\L, ^KWVfW^^ \\u\\Hs-,+K\\Vf\\H.., . 

As in from this and the usual nonlinear estimates (|2.2|) and (|2.4j) . it 
can be verified that there exists a generic function C (depending only on 
parameters that are considered fixed) such that, 

||/i,hf|Ihi, ^ C{\\{9,£p,£v)\\^s+i){l + \\e{dt + V ■ V)p\\hs + k ||0||^.+2}, 

||/2,Hf|Ihi^ ^ C{\\{e,£p,£v)\\j^s+i){l + \\£{dt + V- V)v\\j^s + fl \\£v\\ , 

II/3,hf|Ih; ^ C{\\{9,ep,ev)\\^s+i){l + \\{dt + v- V)0||^. + k 
Set ip = {6,£p,£v). The key point is that 
m + v-V)i;\\Hs 

^C{m^s+i){l + \\{iyVp,udivv,efiV\,KV^e)\\^^ + \\Q\\j,s}- 

This estimate differs from the one that appears in Lemma 5.14 in PP in that 
the right-hand side does not involve v itself but only its derivatives. Yet, as 
the reader can verify, the same proof applies since we do not estimate dt^ 
but instead dtip + v ■ V^p. 

Estimate for fq and f^. By using the elementary estimate 
II (/ - Jsi,)u\\jja+i < ev ||-u||j;^^+i , 

we find that 

\ \\P{xi{4>)Q)\\hi^ + mxmQ)\\Hi ^ \\xi{ct>)Q\\Ht^^ + IIx3(0)QIIh;+^ . 

The tame estimates (|2.2|) and (|2.4|) (see also Lemma 5.5 and 5.6 in ^) imply 

\\x^{m\Hi+^ < (1 + - x^(o)||^;+0 IIQII^.+i < cm\Hi+^) \\Q\\m+' 

so that ||/i,q||2,^(j^ij + \\f'i,Q\\Lf{Hl) ^ ^WQWifim+^y The technique for 
estimating is similar; we find that ||/i,xIIl°°(hi ) ^ ^• 

By definition of |H|;fs(j'), the previous estimates imply that there exists 
a constant C depending only on J7 such that 

r \\iF,,F2)\\H. + m\Hi dt ^Vt( r \\iF,,F2)\\l^ +\\F4l,dty^' 
Jo ^Jo " / 
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From here we can parallel the rest of the argument of Section 5 in [J , to 
prove that \\{Pp, Pv, Pe)\\^^T ^ Cq exp(^/rC) + VTC \\Q\\^^^j^s+i^ where 

the norm IHI^^- is as defined in p.Sp . Since < \\{Pp, Pv, P9)\\^j,, this 
completes the proof. □ 

Step 5: Low frequency estimates. The following step is to estimate the 
low frequency part of the fast components: 

O^F ._ ||div Jei,f ll^oc(ji^s-i) + ||div j£,,7;||^2^(j:^,-) 
+ \\'^JsuP\\l^^{H=-'^) + ^ \\'^'^suP\\l2,(h=) ■ 

Proposition 3.13. Let d ^ 1. There exist constants Cq depending only 
on Mq, C depending only on ft and C depending only on il + T,, such that 

(3.9) ^ Coe(^+^)^ + VtC. 

By contrast with the high frequency regime, the estimate (|3.9|) cannot be 
obtained from the estimates by an elementary argument using differenti- 
ation of the equations (see 12U\ ) . To overcome this problem, we first give 
estimates for the time derivatives, and next we use the special structure of 
the equations to estimate the spatial derivatives. 

For the case of greatest physical interest (d = 3), the proof given in ^ 
applies with only minor changes. Indeed, as alluded to in Remark 13. 9[ it 
suffices to check that all the estimates involving ||w||j|^s remain valid with 
replaced by ||f H^^oo + ||Vt;||j:^s-i. Yet, if d ^ 2, because of the lack of 
estimates for the velocity, we cannot use the time derivatives. For this 
problem, we use an idea introduced by Secchi in |^. Namely, we replace dt 
by the convective derivative 

D^:=dt+v V. 

For the reader convenience, we indicate how to adapt the three main calculus 
inequalities in pi when dt is replaced by D^. 

First, to localize in the low frequency region we use the following commu- 
tator estimate. The think of interest is the gain of an extra factor e. 

Lemma 3.14. Given s > 1 + d/2, there exists a constant K such that for 
all £ G [0, 1], all V G [0, 2], all T > 0, all m gN such that 1 ^ m ^ s and all 
f,u andv in C°°([0, T]; F°°(D)), 

\\[f,Jeu{eD 

m—1 

^ Ke{\\f\\H. + ^ \\A~,'{eD,YD,f\\Hs-^-,] 

m—1 

X {\\A~l^{eD,ru\\^s-^ + ||A-J(eZ),)^n||^._i_.}, 
where K% := {I - {evfAY'"^. 
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To apply the previous lemma, we need estimates of the coefficients / 
and Dvf- Since, for System (|1.4j) . the coefficients are functions of the slow 
variable {9,ep,ev), the main estimates are the following. 

Lemma 3.15. Let s > 1 + d/2 be an integer. There exists a function C(-) 
such that, for all a = (e, n, k) £ A, all T > and all smooth solution 
{p,v,9) £ C°°([0,r];i?~(D)) of dLll), ifu£ [(^ + k)/2,2] then the function 
^ defined by 

:= (^tp, Dyip,'Vip) where ijj := {6,ep,ev), 

satisfies 

(3.10) Yl \\K.HeD,YnH^-i-i ^ C(||^||^._i + S), 

(3.11) WKuieD^Ynnt-' ^ C(||vl/||^._, + S) H^'H^. , 

where S is as defined in (|3.3|) . 

Once this is granted, we are in position to estimate the commutator of 
the equations ()1.4() and V := Je^i^Dv 



/i,LF= [9m,V]D,p + gi{^)[v,V] - Vp - -^[Bi{cP,d^),V]e, 
/2,LF= [g2{4>).V]D,v + g2{(l))[v,V] -Vv - ii[B2{cl),d^),V]v, 

fz,i.Y = [9m,'P]D,e + gi{(i))[v,v] ■ve-K[B^{(i),d^),v]e. 

It is found that 

II/i,lf|Ihi +II/i,lf|Ihi +||/i,lf||^i^(i + ||^|IhOC(||^|Ih-i + s). 

Note that ^ is estimated by means of (|3.8|) . 

As in the high frequency regime, we have to estimate source terms of the 
form e~^VF{'^ , Q). The fact that these large source terms cause no difficulty 
comes from the fact that e'^ Jeu{eD^Y F{m ,Q) = J^y{eD^Y-^D^F{^ ,Q) 
together with D^F{^ ,Q) = 0(1) (the norm S introduced in (|3.3j) is the 
requested norm to give this statement a precise meaning). 

With these results in hands, one can estimate j£y{eD^Y{p, v,9) by means 
of Theorem 13.101 Next, we give estimate for div Je^,?; and VJ^uP from the 
estimate of Jsu{£F)yY{p, v, 0) by means of the following induction argument: 

Lemma 3.16. Set H'uHy^^o-^ji-) := ||^||/,oo^g rp.^a-i-^ + v ||^||/^2^q xw^) ■ 
Let U := {p,v,6) solve 

' 9ii<P){dtP + v-Vp) + e-^ div V - Ke-^xii<P) dw{k{e)Ve) = /i, 
(3.12) <! g2i<l)){dtv + vVv) + e-^Vp- fiB2i<l),d,)v = f2, 
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If support of the Fourier transform ofU is included in the ball {|^| ^ 

then there exist constant Cq depending only on Mq and C depending only 

on i7 such that, for all a G [1, s], 



(3.13) 



k:Z+\t) 

+ eC ||(/i, /2)||,cj(r) + Wf^WiKH") ' 
where C := Coe^^+^^^ . 

Step 6: estimates for the slow components. To complete the proof 
of (|3.2|) . it remains to estimate curlu and 9. Yet, this is not straightforward. 
Following Metivier and Schochet l20!, we begin by estimating curl(7f) for 
some appropriate positive weight 7 = T{6,ep). 

Lemma 3.17. Let d ^ 1. There exist constants Cq depending only on Mq 
and C depending only on fl, and there exists a function T G C°°(M^) such 
that, with 7 = r{6,£p), there holds 

||curl(7w)|| ioo(^.-i) + ^ ||curl(7w)||i2^(^.) ^ Coe^*^ + VTC \\Q\\ L^(Ht+^) ■ 

Lemma 3.18. Let d ^ 1. There exist constants Cq depending only on Mq 
and C depending only on Vt, such that 

||Jei.6'|lioo(^.+i) + ^/ll\\Jeu9\\L2^^Ht+^') ^ CQe"^^ + \/TC \\Q\\ ^^s+i^ . 

The proofs of Lemma l3. 171 and 13. 181 follow from a close inspection of the 
proofs of Lemma 6.25 and 6.26 in 1 . We just mention that this is where 
we use the function F of Assumption (1^1) in ^1.21 (7 is related to the fluid 
entropy) . 

Lemma 3.19. Assume d ^ 3. There exist constants Cq depending only on 
Mq and C depending only on such that, with 70 = T{9q,£Pq) where T is 
as above, there holds 

||curl(7o?;)||^op(^.-i)+V^||curl(7o?;)||^2(j:^.) ^ Cqc'^'^^+VtC \\Q\\^^^jjs+i^ . 

Proof. Set 7 := 7 — 7o- By Lemma [3. 171 all we need to prove is that 
(3.14) ^ 

||curl(7?;)||^op(^.-i) + ||curl(7?;)||^2^(^,) ^ VtC + VtC \\Q\\^^^jjs+iy 
To do so, we claim that 7 is small for small times: 
(3-15) II7|Il-(//.) + II7|Il2,(//»+i) ^ VtC + VtC ||Q|lioc(^s+i) . 
Let us assume this and continue the proof. 
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We have to estimate curl(7w) = 7 curl u + (V7) x v. By combining the 
Cauchy-Schwarz estimate with the usual product rule (|2..3|) and the product 
rule (|2.7|) . we find that 

||7curlf ^ \\i\\l^(h=-i) l|V7;||iop(^.-i) , 

The claim ()3.15() then yields the desired bound H3.14() . 

We now have to prove the claim p.lSf) . We first note that 

^^II7|Il2(h»+i) ^ i^^^II7|Il-(//.+i) 

^ uVfC{\\{e,ep)\\L^^L^^){l + \\{e,ep)\\L^ 

^VTC{\\{9,ep)\\^^^^s+.^)^VTC. 

To prove the second half of (|,3.15|) . we verify that, directly from the defini- 
tions, 7 satisfies an equation of the form dt^ + v ■ V7 = / with / bounded in 
L^(0, T; H^{R'^)) by a constant depending only on Q + \\Q\\l°°(^h;,+^)- Then, 
we apply the above mentioned estimate for hyperbolic equations: 

(3.16) imL^iHs) < e^"" WmWns + £ ||/||^. dt, 

where V = K \\Vv\\jjs-i dt with K = K{s,d). Since 7(0) = 0, by apply- 
ing the Cauchy-Schwarz inequality, it is found that the L'^{H^) norm of 7 
is estimated by \fTe^^ \\f\\L\(H^)^ thereby obtaining the claim. □ 

Step 7: closed set of estimates. To complete the proof of Proposi- 
tion [321 it remains to check that we have proved a closed set of estimates. 
The obvious estimate ||w||j^<t ^ || J£,yw||j:^CT + ||(/ — Jev)u\^„ implies that 

II (Vp, div v)\\ L^^H^'^) + II (Vp, div ?;) || ^2 (j:^^) < I^lf + S^hf, 

and, similarly, ||^||^oc(j:^s+i-) + II^^IIl2 (_f/''+^) estimated by 

II Jei.6'|liop(^j+i) + ll«^ei^V6'||^2 (^2+1) + J^HF- 

The estimate ||eu||j|^<T+i < ||e^i||j;^2 + || Vti||j:^a-i -|- ||(/ — Jeu)'u\\^cT+i yields 

+ ||(Vp,Vt;)||^oo(^.-i) -F^||Vu||^2^(^.) + Ohf- 

On the other hand. Corollary 12.51 implies that, if d 7^ 2, there exists a 
constant Co depending only on Mq such that 

\\^v\\l^{h--^) + \/^ II 

^ Co II (div v,curl(7ov)) 11^00(^^.-1) CqvT^ || (divt;, curl(7oz;))||^2 (j:^,) . 
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By using the estimate (|:-{.8|) . one can verify that the term \\{ep, £v)\\^oo(^^2^ (in 
the left-hand side of H3.17|) ) can be estimated as in the proof of Lemma b.lll 
Therefore, according to Propositions l3.12fiXT^ and Lemma l3.18fiXT^ we 
have proved that, if d / 2, then O < C where C = Coe(^+^)'^ + VTC 
for some constants Cq, C and C depending only on Mq, Q and + S, 
respectively. 

This concludes the proof of Proposition 13.61 and hence Theorem 13.31 
4. Uniform estimates in the Sobolev spaces 

With regards to the low Mach number limit problem, we mention that the 
convergence results ^ proved in apply for general systems (not only for 
perfect gases). To avoid repetition, we only mention that one can rigorously 
justify the low Mach number limit for general initial data provided that one 
can prove that the solutions are uniformly bounded in Sobolev spaces (see 
Proposition 8.2 in pQ). The problem presents itself: Theorem 11.21 onlv gives 
uniform estimates for the derivatives of p and v. In this section, we give 
uniform bounds in Sobolev norms. 

Theorem 4.1. Let d ^ I and N 9 s > 1 + d/2. Assume that Q = 0. 
Also, assume that either xi = Xi(^)P) ^-s independent of or that d ^ 3. 
Then, for all Mq > 0, there exists T > and M > such that, for all 
a = (e,/i, k) G a and all initial data {po,vo,9o) G H^^^{W^) satisfying 

\\{po,vo,Oo)\\fjs+i ^ Mo, 

the Cauchy problem for (11.4(1 has a unique classical solution (p, v, 9) in 
C^{[Q,T]-H'+^{W^)) such that 

sup \\{p{t),v{t),e{t))\\Hs^M. 

iG[0,T] 

The first half of this result is proved in 01. Indeed, the assumption that 
Xi(i?, p) does not depend on -d is satisfied by perfect gases. So we concentrate 
on the second half (d ^ 3). In view of Theorem I3.3| it remains only to prove 
a posteriori uniform estimates. More precisely, the proof of Theorem 14.11 
reduces to establishing the following result. 

Lemma 4.2. Let d ^ 3. Consider a family of solutions {p"',v"', O"") of (|1.4() 
(for some source terms Q"") uniformly hounded in the sense of the conclusion 
of Theorem VJ.iA 

(4.1) sup||(p^^;^nlU^(T) <+°0' 

for some s > 1 + d/2 and fixed T > 0. Assume further that the source terms 
Q"^ are uniformly bounded in C^{[0,T] - L^ n L^(M'^)) and that the initial 
data {p"' (0) , v"" (0)) are uniformly bounded in L'^{W^). Then the solutions 
(^pO-^ya^Qa^ ore uniformly bounded in C°{[0,T]; L^{R'^)). 

"'^These results are strongly based on a Theorem of Metivier and Schochet [2nj about 
the decay to zero of the local energy for a class of wave operators with variable coefficients. 
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Remark 4.3. We allow Q"' for application to the combustion equations. To 
clarify matters, we note that one can replace (|4.1j) by 

sup sup ||(Vp'^(t),V?;'^(t))||^. + ||r(t)||^.+i <+oo, 

aeAte[0,T] 

for some s > 2 + d/2. 
Proof. For this proof, we set 
R := sup{||(p^^;^r)||;,.(^) + IKp'^W, ^;'^(0))||^. + \m\cH[o,T];L^nL^)} ^ 

and we denote by C{R) various constants depending only on R. 

The strategy of the proof consists of transforming the system ()1.4() so as to 
obtain estimates uniform in e by a simple integration by parts argument. 

To do that we claim that there exist V £ C^([0, T]; L2(M'^)) satisfying 
the following properties: 

(4.2) sup||(p^^;'^)||^^(^.)^sup||^'^||^«,(^.) + C(i^), 

(4.3) sup||[/'^(0)||^. ^C(i?), 

aeA 

and 17°' solves a system having the form 

(4.4) E^dtU" + e-^Sid,^)^ = F", 

where S{dx) is skew-symmetric, the matrices E"" = E°'{t,x) are positive 
definite and one has the uniform bounds 

(4.5) sup||9i£^''||^cx,([o,2^]xMd) + ||(-E'')^^||^^^jQ^yj^jgd^ + ||F"||^i^ ^ C{R). 

Before we prove the claim, let us prove that it implies Lemma l4.2l To see 
this, we combine two basic ingredients: 

^( S'^t/" = 5(a^)[/^ [/'^ ) + ( F^ [/'^ ) + ( (atS")C/" , ) 

^ \\F''\\l2+C{R)\\U''\\l,, 

and WU'^Wli < \\{E'')-m'^l{E''U'',U''). Hence, by and the 

Gronwall's Lemma implies that ||f^"||j;^cx)(j;^2) ^ C{R). The estimate (|4.2p 
thus implies the desired result. 

To prove the claim, we set C/" := {p"',v"' — V"')'^ where 

V ■= KXi(0")A:(0")Vr + VA-1(-kVxi(<^") • A;(r)Vr + xi('/''')Q''). 

The fact that V"" is well defined follows from Proposition l2.11 We next verify 
that satisfies with 

j^a_(gim \ ^.^^_/^0 div 



/ -gi{r)v''-Vp'' 
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By dSini), (El) and to prove that the bounds and hold, it 

suffices to prove that ^ C{R). Yet, this is nothing new. Indeed, 

we first observe that, directly from the equations. 

On the other hand, the product rule H2.7|) implies that Hu"^ • \/(j)°'\\^3-i is 
estimated by ||Vt'"||j^s-i ||</'"||//s ^ C{R). This completes the proof. □ 

Remark 4.4. For our purposes, one of the main differences between 
and M is the following. For all / G Cg°(M^), Proposition 12.51 implies that 
there exists a vector field u G H°^{M.^) such that dWu = f. In sharp 
contrast, the mean value of the divergence of a smooth vector field u G 
H°°(R) is zero. This implies that Lemma 14.21 is false with d = 1. 

The following result contains an analysis of the easy case where initially 
9o = 0{e). This regime is interesting for the incompressible limit (see 

Proposition 4.5. Let d ^ 1 and R 3 s > 1 + d/2. For all Mq > 0, 
there exists T > and M > such that for all a £ A and all initial data 
iPo,vo,9o) e if*(M^) satisfying 

(4.6) ll(Po,^^o)ll//»+e"'ll^oll//= ^Mo, 

the Cauchy problem for (|1.4() has a unique classical solution {p, v, 9) in 
C°([0,r];i?*(R"')) such that 

(4.7) sup mt),v{t))\\H^+e-^ \\0{t)\\Hs ^ M. 
te[o,T] 

Proof. The proof of this result is based on the change of unknown (p, v,9) 
{G{9,ep),v,9) where G is as given by Assumption (lEJ in ^1.21 By setting 
p = G{9, ep) it is found that (p, f , 9) satisfies ()1.4() if and only if 

' X'iidtP + v ■ Vp) + {x3 - xOdivw = 0, 

(4.8) I g2{dtv + vVv) + e-^jiV9 + e-^-f2^p - fiB2V = 0, 

, 93{dt9 + v ■V9) + diYv - kx3 diY{kV9) = 0, 

where 71 = (xids) / {X39i) 72 = l/ffi- Notice that Assumption (HU 
implies that the coefficients gi, 'ji, X3 and X3 ~ Xi are positive. 

The key point is that the assumption ()4.6|) allows us to symmetrize the 
equations by setting u := {p,v,9), where 

p:=e~^p, 9:=e-^9. 

The fact that this change of unknowns is singular in e causes no difficulty. 
Indeed, directly from the assumption (|4.6|) . we have ||^(0)||j:^s ^ Mq. On 
the other hand, the assumption G(0, 0) = implies that there is a function 
Cg such that 110(^)11 j^^,, ^ CcdluHioo) \\u\\h'^ for all u e with a > d/2. 
Therefore, we have 

\\p\\Hs=e-'\\G{e,ep)\\Hs^e-'CG{\mep)h^)\\{9,ep)\\Hs 

^CG{\\{e,ep)\\^^me,p)\\Hs, 
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hence, ||/3(0)||j:^s ^ Co for some constant depending only on Mq. 

Because p) ^ {'d,G{'d,p)) is a C°° diffeomorphism with G(0,0) = 
0, one can write ep = P{9,G{9,epj) = P{9,p), for some C°° function P 
vanishing at the origin. Therefore one can see the coefficients {gi, Xi^ li---) 
as functions of {0,p). Hence, with u = {p,v,0) as above, one can rewrite 
System (|4.8j) under the form 

(4.10) Ao{eu)dtu+ ^ Aj{u,eu)djU+- ^ Sj{eu)djU—B{eu,dx)u = 0, 

where the matrices Sj, Aj are symmetric (with Aq positive definite) and the 
viscous perturbation B{eu,dx) is as in 1)4. 8() . Furthermore, one can always 
assume that the matrices Sj have constant coefficients. 

Since the matrix Aq multiplying the time derivative depends only on the 
unknown through eu, and since the initial data u{0) are uniformly bounded 
in H^, the proof of the uniform existence Theorem of [T^ applies. By 
that proof, we conclude that the solutions of (|4.inj) exist and are uniformly 
bounded for a time T independent of e. Once this is granted, it remains 
to verify that the solutions (p, v,9) of System (|1.4|) exist and are uniformly 
bounded in the sense of (|4.7() . To see this, as for p in (|4.9() . we note that 

\\p\\hs = \\P{0,p)\\h^ 

^e~'Cp{\\{9,p)h^)\\{9,p)\\Hs=Cp{\\{9,p)\\L^me,p)\\Hs 

^C{\\{9,p)U^), 

so that ||(p, f )||j|^s + II^IIhs ^ CdluHji^s). This completes the proof. □ 

Remark 4.6. Consider the Euler equations (/x = = k and e = 1). By 
a standard re-scaling. Proposition 14.51 iust says that the classical solutions 
with small initial data of size 5 exist for a time of order of 1/5. Following 
the approach initiated by Alinhac in J^, several much more precise results 
have been obtained. In particular, the interested reader is referred to the 
recent advance of Godin (for the 3D nonisentropic Euler equations). 

5. Spatially periodic solutions 

In this section, we consider the case where x belongs to the torus T'^. 

Theorem 5.1. Let d ^ 1 and N 9 s > 1 + d/2. For all source term 
Q e C°°(M X T"^) and for all Mq > 0, there exist T > and M > such 
that, for all a £ A and all initial data {po,V(),9q) G H^^^{T'^) satisfying 

\\{po,vo)\\h^ + \\{0o,£Po,£Vo)\\hs+i ^ Mo, 

the Cauchy problem for (|1.4() has a unique classical solution {p, v, 9) in 
C°([0,r];i7*+^(T^)) such that 

sup ||Vp(t)||^._i + Wvmn^ + \mt),epm\H^ < M. 

t€[0,T] 
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The proof follows from two observations: first, the results proved in 
Steps 1-6 in section |21 apply mutatis mutandis in the periodic case; and 
second, as proved below, the periodic case is easier in that one can prove 
uniform estimates for the velocity. This in turn implies that (as in ^ [20] ) 
one can directly prove a closed set of estimates by means of the estimate: 

||f llHs(Td) ^ C \\divv\\jjs-i(^Td) + C \\cml{jv)\\ fjs~i(jd-^ + C \\v\\i^2(jd) , 

for some constant C depending only on ||log7||j|^sC]pd) (compare with (|2. 151) 1. 
Let us concentrate on the main new qualitative property: 

Lemma 5.2. Let 1. Consider a family of solutions {p"" ,v'^ ,6"") of (|1.4j) 
(for some source terms Q"") such that 

sup||(p^^;^r)||;^.(^) < +CX), 

for some s > 1 + d/2 and fixed T > 0. // is uniformly bounded in 
Ci([0,r];L2(T'^)) and (p"(0), w'^(O)) is uniformly hounded in ^^(T^), then 
v"" is uniformly hounded in C'^([0, T]; L^(T°')). 

Proof. The main new technical ingredient is, as used by Schochet in |23j . an 
appropriate ansatz for the pressure. 

Again, the proof makes use of the Fourier multiplier VA~^. Note, that 
VA-i is bounded from L|(T'^) to H^{T'^) where L^(T'^) consists of these 
functions u G L^(T'^) such that (u) := /j.^ u{x) dx = 0. 

Set 

and introduce the functions V"" = V"'{t, x) and P"^ = P°-{t) by 

^"•=7-^Sv and :=VA-i(F"-5i(</)'^)P"), 

so that 

F" = 5i(0'')P'' + divy^ 
This allows us to rewrite the first equation in 1)1.4(1 as 

gi{<l,''){dtp'' + ■ Vp'') + e-^iviv'^^ - = gimP"- 
Therefore, by introducing 

[/« := (g'^,^;"-^'^)^ with q''{t,x)=p''it,x)-P''{t), 
we are back in the situation of Lemma 14.21 U"" satisfies 
(5.1) E'^idtU" + v" ■ VU") + £-^3(8^)^ = F^ 

where S{dx) is skew-symmetric, the matrices E"" are positive definite and 

\\{E°- , 8tE°- + v"" • VF")||^oo([o,T]xKd) + ||(^")~^ILi([o,T]xR'') + ll-^^llii^CLa) 

is uniformly bounded. 
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As before, the proof proceeds by multiplying by U"' and integrating on 
[0,r] X T'^. We find that dtiE^U" , U" ) is given by 

{{{dt + v''- V)E'') [/^ C/'^ ) + ( ^'^(div v'')U'' + 2( F^ [/'^ ), 

and hence conclude that C/" is uniformly bounded in C^{[0,T]; L'^iT'^)). 
Since V"- is uniformly bounded in C'^([0, T]; L^(T'^)), this yields the desired 
result. □ 

Remark 5.3. In the periodic case, as shown by Metivier and Schochet \21\ 
\22\ as well as Bresch, Desjardins, Grenier and Lin jHI, the study of the 
behavior of the solutions when e ^ involved many additional phenomena. 

6. Low Mach number combustion 

The system is relevant whenever all nuclear or chemical reactions are 
frozen, which is the case in many treatments of fluid mechanics. By contrast, 
for the combustion, one has to replace the energy evolution equation by 

dtipe) + div(p?;e) + Pdivt; = Kdiv(A;Vr) + F{Y), 

with Y := (Yi, . . . ,Yl) where the Y^s denote the relative concentrations of 
nuclear or chemical species. The new unknown Yi satisfies : 

(6.1) dtipYe) + div{pvYe) = Xdiv{DiVYe) + poje{t,x), 

where ivg is a given source term, > and A measures the importance of 
diffusion processes. 

Many results have been obtained for the reactive gas equations (see 
and the references therein). Yet, the previous studies do not include the 
dimensionless numbers. Here we consider the system: 

( a{dtP + V ■ VP) + dwv = K(3dW{kVT) + Fi{Y,T, P), 



(6.2) 



VP 

p{dtv + V ■ Vv) + = p[2dw{CDv) + V(r/divi;)), 



-f{dtT + V ■ VT) + divt; = K6diY{kVT) + F3{Y,T,P), 

^ p{dtY + V ■ VY) = Adiv(DVy), 

where a, /3, 7 and 5 are given functions of {Y-, T, P). 

As explained in the introduction, it is convenient to introduce {p, 9, y) by 
p = Pe'^P, T = T/, Y = Yey, where {P,T,Y) G [0,+oof+^. For smooth 
solutions, {p,v,0,y) satisfies a system of the form: 
(6.3) 

gii'P)idtp + i; • Vp) + ^ divi- = ■^xi(^) dW{k{e)Ve) + ^Qi($), 

g2i^)idtv + vVv) + ^Vp = px2i^){div{C{9)Dv) + V{t]{9) divt;)), 
g3i^)idt9 + vV9) + divv = Kxsi'^) div(A:(^)V^) + Q3($), 
, gmidty + v-Vy) = Xx4(.^) dw{D{9)Vy), 
where <I> = {y, 9, ep). 
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Assumption 6.1. Denote by (y, i?, p) G M.^ the place holder of the unknown 
{y,6,ep). Parallel to Assumption (lEJ in ^1-21 we suppose that gi and Xi 
{i = 1,2,3) are positive functions of (y,'(?,p) G M^, xi < Xs and 
there exist two functions F and G such that (y, "!?, p) i-^ {F{y,-d, p), p) and 
(y, 1?, p) (y, ??, G('(?, p)) are C°° diffeomorphisms from onto R-^, F and 
G vanish at the origin, and 

dF dF dG dG 

ov op ov op 

Moreover, Qi and are C°° functions of (y, "i?, p) vanishing at the origin. 
Introduce 

B := { (e, fi, K, A) G (0, 1] x [0, 1] x [0, 1] x [0, 2] | A ^ ^/J^^7^}■ 

Definition 6.2. Let T > 0, s G M, 6 = (e, /x, k. A) £ B and set a := (e, fi, k) . 

The space Z^{T) consists of these {p,v,6,y) G G°([0, T]; i?"(]R^)) such that 

{p,v,e)eX!{T), uyeG'{[0,T];H^+\R^)), Xy e L\0,T; HI+\R^)), 

where v := ^/J^+l^ and '^^(T') is as defined in Definition The space 

2^{T) is given the norm 

\\{p,v,e,y)\\zs^T) ■■= ll(p>^>^)IU^{T) + + ^llyllL2(H,?+2) • 

Having proved estimates for the solutions of System (|1.4|) with precised 
estimates in terms of the norm S of the source term Q (see we are 

now in position to assert that: 

Theorem 6.3. Assume that d / 2. Given Mq > and N 9 s > 1 + d/2, 
there exist T > and M > 0, such that for all b £ B and all initial data 
{po,vo,9o,yo) G H'+^iR"^) satisfying 

\\{Vpo,Vvo)\\hs^i + \\{yo,Oo,epo,evo)\\j^s+i ^ Mq, 

the Cauchy problem for ()1.4() has a unique classical solution (p, v, 6, y) in 
the ball B{ZI{T);M). 

Remark 6.4. For the case of greatest physical interest {d = 3), Theorem l6.3l 
has two corollaries. As alluded to in Section it allows us to rigorously 
justify, at least in the whole space case, the computations given by Majda 
in JH]- By the way, this proves the well posedness of the Cauchy problem 
for the zero Mach number combustion in the whole space (this was known 
only in the periodic case J^). Moreover, note that the solutions given 
by Theorem 16.31 satisfy uniform estimates recovering in the limit e — > 
those obtained by Embid for the limit system. Finally, we mention that the 
previous analysis seems to apply with Qj($) replaced by Xi{^)Q{^, Vy, V^y) 
for some smooth function Q, yet we will not address this issue. 
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Appendix A. General equations of state 

Recall that, in order to study the full Navier-Stokes equations we 
choose to work with the unknown (P,v,T). In order to close this sys- 
tem, we must relate (p, e) to (P, T) by means of two equations of state: 
p = p{P, T) and e = e(P, T). The purpose of this section is to show that As- 
sumption (H2I) in ^1-21 is satisfied under general assumptions on the partial 
derivatives of p and e with respect to P and T. 



A.l. Computation of the coefficients. We begin by expressing the co- 
efficients Qi and Xi) which appear in H1.4|) . in terms of the partial derivatives 
of p and e with respect to P and T. To do that it is convenient to introduce 
the entropy. Here is where the first identity in H1.5|) enters. 

Assumption A.l. The functions p and e are C°° functions of {P,T) G 
(0,-|-oo)^, satisfying 

pdp_ dp_ ^ ode_ 
dP dT ^ dP 

Introduce the 1-form uo defined by Tuj := de-|-Pd(l//9), where we started 
using the notation d / = {df /dT) dT + {df /dP) dP. Assumption E3] im- 
plies that da; = 0. Hence, the Poincare's Lemma implies that there exists 
a function S = S{P,T), defined on (0,+cxo)^, satisfying the second 
principle of thermodynamics: 

(A.l) TdS = de + Pd{l/p). 

By combining the evolution equations for p and e with HA.1|) written in the 
form pTdS = pde — (p/p) dp, we get an evolution equation for S, so that 

+ " • ^) {s) = (T (ptVO Cdiv(wr) + q) ■ 

On the other hand, one has 



Equating both right hand sides and inverting the matrix J, we obtain 
(A.2) 



{dtP + V ■ VP) +adivv- Kb div{kVT) = bQ, 
{dtT + V ■ VT) + cd\Yv- Kdd\\{kVT) = dQ, 



where 



p(dS/dT) , dp/dT p{dS/dP) , dp/dP 

= — — 7-—) C = ; — T 1 d — 



det(J) prdet(J) det(J) prdet(J) 
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To express the coefficients gi and Xi i^i terms of physically relevant quan- 
tities, we need some more notations. We introduce 

^ • ^ 95 ^ (g5/9r)(9;./gP) - {dS/dP){dS/dT) 

^= ^ar' ■= ^ 

The functions Kt, Kp, Cy and Cp are known as the coefficient of isothermal 
compressibility, the coefficient of thermal expansion and the specific heats 
at constant volume and pressure, respectively (see Section 2 in J2])- The 
function TZ generalizes the usual gas constant: for perfect gases one can 
check that IZ = R. 

We now have to convert System IA.2I into equations for the fluctuations p 
and 9 as defined by (|1.2|) . Performing a little algebra we find that 



^^{d,p + . . Vp) + 1 div. - div(fcrV^) = i^^Q, 

Cp £ £ pCp £ pCp 

p{dtv + v- Vv) + ^PVp = /u(2div(C-Dw) +V{r]dwv)), 
^ pCvT{dte + vV9)+ TZpT div v - Kdiv{kTVe) = Q. 
Hence, {p, v, 9) satisfies (|1.4|) with 

* KtCvP , p , Cv , Kp , 1 , 1 
(A.4) 5i = — , = d' 53 = Xi = -7^-' X2 = d' Xs 



Cp p Ti pCp p UpT 

where we used the following notation: for all /: (0, +oo)^ M, 
(A.5) r(^,p):=/(^e^PeP). 

A. 2. Properties of the coefficients. 

Assumption A. 2. The functions p and e are C°° functions of (P, T) £ 
(0, +oo)^ such that, p > and 

, . , dp dp , de dp de dp 

Remark A. 3. This assumption is satisfied by general equations of state. 
Indeed, ()A.6|) just means that the coefficients Kt, Kp and Cy are positive. 



The following result prove that Assumptions I A. ll and [^31 implv that our 
main structural assumption is satisfied. 

Proposition A.4. If Assumvtions \A. 1\ and \A.Sl are satisfied, then xi < Xs 
and Qi, Xi (^ = 1,2,3^ are C°° positive functions. 

Proof. In view of (|A.4|) , the proof reduces to establishing that 

C 

0<Kt, 0<Kp, 0<Cv <Cp and < 7^ < — 

1 Kp 



LOW MACH NUMBER FLOWS 27 

The first two inequalities follow from the definitions of Kt and Kp. To prove 
the last two, we first establish the Maxwell's identity dS/dP = p~'^{dp/dT). 
To see this, by ()A.1|) . we compute 

|^dTAdP = d(Td5)=d{de + Pd(i)} = --i^dPAdr. 

Since dp/dT < 0, the Maxwell's identity implies that dS/dP < 0. By 
combining this inequality with dp/dP > 0, we find TZ > 0. Also, the 
identity dS/dP = p^'^{dp/dT) imphes that 

Cp _ {dS/dT){dp/dP) 

Cv ~ {dS/dT){dp/dP) - p-^{dp/dTf' 

which proves Cy < Cp. 

In view of (lA.ll) . the assumption — — -— - > -rr-—— is equivalent to 
^ ^' ^ dTdP dPdT ^ 

dS dp dS dp 

dfdP ^ 'dPdf' 

This inequality has two consequences. Firstly, it implies that Cy > 0. 
Secondly, it yields 

TKpTZ _ {dS/dP){dp/dT) 
Cp ~ {dS/dT){dp/dP) ^ ■ 
This concludes the proof. □ 

We now discuss the physical meaning of the functions F and G introduced 
in ^1.21 These are compatibility conditions between the singular terms and 
the viscous terms. To see this, suppose (p, v, 6) is a smooth solution of H1.4() 
and let ^' = ^'(i?, p) G C°^(M2). Then ^ := '^{0,£p) satisfies 

gm {dtip + V • VV') 

+ (51^ + 53^) divt; = k(5iX3^ + 53Xi^j(div(Me)Ve) 

= :ri(*) =:r2{*) 

where the coefficients gi, Xi-, d'^/d'd and d"^ /dp are evaluated at {6,ep). We 
next show that for appropriate function ^ one can impose 

(A.7) [ri(^') = and r2(^) > 0] or [ri(^) > and r2(^) = 0]. 

Proposition A. 5. Assume that Assumvtions \A.l\ and lX'^ are satisfied and 
use the notation ()A.5|) . The functions S* and p* satisfy 

, dS* dS* dp* dp* 

Remark A. 6. The fact that ^' = 5* (or ^' = p*) satisfies the first (respec- 
tively second) set of conditions in ()A.7|) now follows from xi < X3- 
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Proof. By (|A.4|) and the definitions given in ()A.3|) . one has 

(A 9) g!_ P{dS/dP) 

^ ' gl T{dS/dT) 

By definition (H^, df*/d^ = [T{df/dT)Y and df*/dp = [P{df/dP)Y. 
This proves that S* satisfies the first identity in (|A.8|) . Next, we compute 

xl _ {dp/dT){ds/dP) 

Xl {dp/dP){dS/dT) 

By (jXH), this yields x\9lP{dp/dP) = -xl9*iT{dp/ dT). Which proves that 
p* satisfies the second identity in (|A.8() . □ 

Remark A. 7. Assumption (Id) in ^1.21 requires, in addition, that F = S* 
and G = p* define bijections. This means nothing but the fact that the 
thermodynamic state is completely determined by (P, T), or (P, S) or (p, T). 

The following result contains an example of equation of state such that 
Xl depends on -d. 

Proposition A. 8. Assume that the gas obeys Mariotte's law: P = RpT , 
for some positive constant R, and e = e{T) satisfies Cy := de/dT > 0. 
Then, Assum'ptions \A.l\ and VK^ are satisfied. Moreover, 

xl = R/{{Cv{T) + R)P), 

so that Xi{'^^ p) 'is independent of 'd if and only if Cy is constant. 
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